Symmetry breaking effects upon bipartite and multipartite entanglement in the XY 

model 
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We analyze the bipartite and multipartite entanglement for the ground state of the one- 
dimensional XY model in a transverse magnetic field in the thermodynamical limit. We explic- 
itly take into account the spontaneous symmetry breaking in order to explore the relation between 
entanglement and quantum phase transitions. As a result we show that while both bipartite and 
multipartite entanglement can be enhanced by spontaneous symmetry breaking deep into the ferro- 
magnetic phase, only the latter is affected by it in the vicinity of the critical point. This result adds 
to the evidence that multipartite, and not bipartite, entanglement is the fundamental indicator of 
long range correlations in quantum phase transitions. 
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I. INTRODUCTION 

One of the most remarkable features of many phase 
transitions is the occurrence of spontaneous symmetry 
breaking, in which the symmetry of the Hamiltonian is 
not realized in the system state. When this occurs a 
macroscopic observable (the order parameter) emerges, 
which is required for a unique specification of the micro- 
scopic state Those quantum phase transitions (QPT) 
occur at zero temperature, and are triggered by the varia- 
tion of a parameter of the system's Hamiltonian 0|. The 
system's eigenenergies then show non-analytical behavior 
which embodies the order of the phase transition. These 
non-analyticities in turn are reflected in several macro- 
scopic observables. 

Lately there has been an increasing interest in describ- 
ing QPTs, not by means of the non-analyticities of the 
spectrum or of the physical observables, but rather by 
the amount of entanglement (bipartite or multipartite) 
present in each of the system's phases. This is moti- 
vated in particular by a general expectation that entan- 
glement can be given as important a status as the energy: 
both quantities can be seen as resources useful for the ac- 
complishment of interesting physical tasks [jj. Moreover, 
long-range correlations often found in strongly correlated 
many-body systems at zero temperature have a purely 
quantum origin and are expected to be inextricable from 
entanglement 0, [![• Several authors have studied the 
role of entanglement in QPTs by considering either bi- 
partite or multipartite entanglement measures (e.g., see 
[6]-[36]), both calculated for spin-1/2 lattice models such 
as the Ising and XY model in a transverse magnetic field 
[37l . [H, [M, S3|- Invariably though many of these de- 
velopments do not take into account the spontaneous 
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symmetry breaking accompanying the QPT, employing 
instead symmetric states in calculations in the ordered 
phase (for an exception see [H, US] and the note at the 
end of the manuscript). This procedure, although com- 
mon and even correct for finite systems, is unrealistic in 
the study of QPTs since it is well known that symmetric 
states (Schrodinger cats) are never realized in the ther- 
modynamic limit due to superselection rules [l[. 

Recently, we have developed a program on the investi- 
gation of multipartite entanglement properties in QPTs 
by propo sing a Generalized Global Entanglement (GGE) 
measure [4ll.l42l|. As the name suggests this is a general- 
ization of the Meyer and Wallach Global Entanglement 
[43 | (here dubbed G(l)) to account for all the possible 
bipartitions of the system state. Besides being able to 
detect any kind of entanglement present in the system 
[12], this measure is also able to signal the location and 
the type (order) of QPTs [H[. Our results developed for 
the infinite 1-D Ising and XY chains indicate that mul- 
tipartite entanglement is maximal at the critical point, 
playing a major role in the QPT process, in contrast to 
bipartite entanglement 0, 0|. 

To correctly consider the spontaneous symmetry 
breaking due to quantum fluctuations is an essential el- 
ement for the success of such an entanglement measure. 
However this point is particularly unclear in the recent 
literature on entanglement in QPTs. Many authors are 
not specific in their choices and have erroneously applied 
a symmetric state in their investigations. In this article 
we extend the entanglement analysis for the XY model 
taking into account the role played by spontaneous sym- 
metry breaking in both bipartite (pairwise) and multi- 
partite entanglement. As we argue in the body of the 
paper, symmetry breaking favors multipartite entangle- 
ment. While bipartite measures (concurrence and nega- 
tivity) are the same irrespective of the state one is em- 
ploying, multipartite entanglement is not. We show that 
both G(l) as well as G(2,n), an auxiliary function defin- 
ing one class of GGE, exhibit completely different be- 
haviors depending on whether the state is symmetric or 
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not. We also review some of the results on bipartite and 
multipartite entanglement in a clearer and more detailed 
fashion. We expect with this work to settle some is- 
sues concerning the relation between entanglement and 
quantum phase transitions, at least in the context of the 
one-dimensional XY model. This paper is structured as 
follows. In Sec. Owe review the essential results in the 
literature regarding bipartite and multipartite entangle- 
ment present in the 1-D Ising and XY models. In Sec. IIIII 
we discuss the XY model in detail and show how the re- 
duced two-spin state is calculated. In Sec.[lV]we analyze 
bipartite entanglement measures (concurrence and nega- 
tivity) by comparing the results obtained for the symmet- 
ric state with the ones obtained for the broken-symmetry 
one. In Sec. |V] we analyze the multipartite entangle- 
ment as given by G(l) and G(2, 1) for the two choices of 
ground states (symmetric or broken-symmetry). Finally, 
in Sec.[Vl]a discussion ends the paper. 



II. ENTANGLEMENT AND QPT IN 1-D ISING 
AND XY SPIN CHAINS 

In this section we outline some of the most relevant 
findings [45| associated with entanglement in the 1-D 
Ising and XY models. 

The first approach we mention, concerning pairwise en- 
tanglement (concurrence) between two spins in the chain, 
was considered in Refs. @, 0]. It was demonstrated 
that the concurrence between nearest neighbors of the 
XY model is maximal not at the critical point but in 
its vicinity. Furthermore, the pairwise entanglement be- 
tween neighbors more than three sites apart vanishes in 
the quantum Ising chain. The authors of Ref. pj also 
showed that the derivative of the nearest neighbors con- 
currence is able to signal the QPT as it diverges at the 
critical point and exhibits finite-size scaling. Thus, it be- 
came clear that the ability to signal a QPT could be a 
general property of good entanglement measures. 

The first work to establish a formal relation between a 
QPT and bipartite entanglement measures was Ref. [§]. 
The authors have demonstrated that, under a set of rea- 
sonable assumptions, a discontinuity in a bipartite entan- 
glement measure (concurrence [§] and negativity pl|) is a 
necessary and sufficient indicator of a first order quantum 
phase transition (1QPT), which is generically character- 
ized by a discontinuity in the first derivative of the ground 
state energy. Furthermore, they have shown that a dis- 
continuity or a divergence in the first derivative of the 
same measure (assuming it is continuous) is a necessary 
and sufficient indicator of a second order QPT (2QPT), 
which is generically characterized by a discontinuity or a 
divergence of the second derivative of the ground state 
energy. Subsequently, it was pointed out [ll| that this 
result was more general and would apply to any entangle- 
ment measure dependent on the reduced density operator 
of two spins. Finally, it was demonstrated in Ref. [T3|, 
using the Density Functional Theory formalism, that any 



entanglement measure can be expressed as a unique func- 
tional of the set of first derivatives of the ground state en- 
ergy. For most of the cases, however, the explicit expres- 
sion of the functional is not known. This result showed 
that any entanglement measure can in principle signal 
a QPT, since it inherits the non- analytical behavior of 
the derivative of the energy. Of course, depending on the 
definition of the entanglement measure used, "accidental" 
cancellations of such divergences / discontinuities may oc- 
cur (see Refs. pi El). 

Another approach to understand 
pairwise entanglement in QPTs based on the study of 
the crossing of energy levels has also been proposed [la ]. 
For the case of a system of indistinguishable particles it 
was proved that, given some provisos, the entanglement 
between one part (A) and the rest (B) is able to signal 
a QPT pj|. However, in this case parts A and B corre- 
spond to modes not particles, in contrast to the former 
mentioned works. We should also note that pairwise en- 
tanglement in small chains (two, three, and four spins) 
for the XY model was previously studied in Ref. [17| . 

The second kind of approach worth mentioning focuses 
on multipartite entanglement (ME). In Refs. P1E3] the 
entropy of entanglement between one part of the chain 
(a block of L spins) and the rest is employed for this pur- 
pose. There the entanglement entropy is defined through 
the von Neumann entropy of one of the reduced parts, a 
valid approach whenever the global state is pure. It was 
shown for some spin- 1/2 models that at the critical point 
(CP) the entanglement entropy increases logarithmically 
with L, whereas it saturates for large L away from the 
CP, a result which had been known from conformal field 
theory (20L [2l|. For the one-dimensional XY model the 
block entanglement was extensively and carefully stud- 
ied in Ref. [22]. Another approach for ME investigation 
(ME) is considered in Ref. [23| through the study of the 
maximal possible overlap between the state studied and 
all possible separable states; the larger this overlap is, 
the less entangled the state. The XY model was analyzed 
in that way and it was shown that the ME is maximal 
around the CP and its derivatives diverge as the CP is 
approached. It was also shown that the ME is zero at the 
second critical point (2CP) where the state is known to be 
separable [11 (see the discussion in Sec. IIV[) . Tripartite 
entanglement, given in terms of the residual tangle [25| . 
was also analyzed for an Ising chain of 3 spins in a trans- 
verse field in Ref. [26]. In this last work it was shown 
that the residual tangle is not maximal around the ex- 
pected CP, which really exists only in the thermodynamic 
limit. With the purpose of studying ME, a new measure 
was defined and analyzed for spin chains in Ref. [23]. It 
was named Localizable Entanglement and defined as the 
maximal amount of entanglement that can be localized 
in two particles, on average, by doing local measurements 
on the rest of the particles. The Localizable Entangle- 
ment was shown to be maximal at the "critical" point for 
a finite Ising chain of 14 spins. This, together with the re- 
sults for the block entanglement, were the first evidences 
that multipartite entanglement could be important in the 
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context of quantum phase transitions. It was also demon- 
strated that connected correlation functions are a lower 
bound for Localizable Entanglement, a remarkable result 
enabling the system to have a finite correlation length but 
infinite entanglement length (see [28| for an example of 
such behavior). We should remark, however, that it has 
been argued that Localizable Entanglement may not be 
an entanglement monotone [29]. 

With the hope that the tools of Quantum Information 
and Computation could help to better understand Quan- 
tum Phase transitions Zannardi and coworkers [3(J have 
proposed and (in their own words) "showed that quantum 
fidelity - the overlap modulus - of two finite-size ground 
states corresponding to neighboring control parameters is 
a good indicator of quantum phase transitions. Indeed, 
the fidelity typically drops abruptly at the critical points, 
as a consequence of the dramatic state transformation in- 
volved in a transition". For the sake of completeness we 
should also mention studies of the t emp oral evolution of 
the entanglement in the XY chain [31| as well as other 
attempts to show that multipartite entanglement is im- 
portant/enhanced in quantum phase transitions, as can 
be found in references 0, [13, M, H SH] . 

We note that none of the employed entanglement mea- 
sures in the above studies are maximal at the CP, with 
the exception of the single site entropy of the Ising model 
[(| in the thermodynamic limit and the Localizable En- 
tanglement [13] of an Ising chain of a few spins. We 
should also mention that in Ref. [3(| the authors have 
studied the loss of entanglement along the renormaliza- 
tion group flow of an XY chain. For this purpose they 
obtained the entanglement between 100 spins and the rest 
of the chain, as a function of the transverse magnetic field 
and the anisotropy, showing that it was indeed maximal 
at the critical point. 

At this point we should mention another interesting 
feature observed in Refs. 0, 0], independently. They 
showed that bipartite entanglement vanishes when the 
distance between the two spins is greater than one lattice 
site. This is quite surprising since long range quantum 
correlations are expected to be present at the CP. It was 
then conjectured that bipartite entanglement at the CP 
would decrease in order for the ME to increase, due to 
entanglement sharing 0. In other words, ME only ap- 
pears at the expense of pairwise entanglement and at the 
CP we should expect a genuine multipartite entangled 
state. 

In Refs. [ill . |42| three of us used the fact that for the 
quantum Ising chain the entanglement between one spin 
and the rest of the chain (given by the averaged linear 
entropy [4(1) is equal to the Global Entanglement (GE), 
G(l), a proposed ME measure introduced by Meyer and 
Wallach in 2002 j43|, in order to show that G(l) is maxi- 
mal at the critical point. Inspired by the GE, we have also 

(n) 

proposed the GGE, E G , where the averages are taken 
over the linear entropy of two, three, and more spins 
(or subsystems). A similar approach was independently 
presented in Ref. [13]. In that construction, we allow the 



spins to be non-contiguous along the chain and not just in 
a continuous block as already considered in Refs. [la.[l9|. 
For example, E G is the entanglement between two spins 
and the rest of the chain averaged over all possible dis- 
tances between the spins. In this context, another quan- 
tity that is also interesting is the average entanglement 
between two spins n sites apart and the rest of the chain 
(without averaging over n): G(2,n). In Refs. [ill. I42I] 
we then analyzed the entanglement between two spins 
n sites apart (G(2,n)) and the rest of the chain show- 
ing that it is maximal at the critical point and increases 
with 71, saturating at the value 0.675 in the limit of large 
n. This result was one of the first indications that mul- 
tipartite entanglement is maximal/enhanced and more 
distributed at the critical point, adding strength to the 
conjecture of T. J. Osborne and M. A. Nielsen [fj|. It 
also suggested that ME is the key ingredient for the ap- 
pearance of the long-range correlations that develop at 
the critical point. In Ref. we have pursued the dis- 
cussion of Ref. [il| further, exploring the advantageous 
features of the Generalized Global Entanglement for an 
operational multipartite entanglement classification and 
quantification, in comparison to the other available mea- 
sures for both finite and infinite collections of two-level 
systems. 

In a more recent work [!i[, we extended the above 
results to the one-dimensional XY model showing that 
G(l) and G(2, n) are maximal at the critical point. In 
a more general context we also showed explicitly that 
G(2,n) is able to signal QPTs, something already ex- 
pected from the results of [ill . ll2T |. Finally, and more 
striking, we demonstrated that for collections of two-level 
systems with symmetry-breaking second-order quantum 
phase transitions, G(2,n) increases exponentially with 
n away from the critical point. This increase is gov- 
erned by a characteristic length, named the entanglement 
length £e, which is half the correlation length £g = £c/2- 
Furthermore, at the critical point G(2, n) increases as a 
power law, implying an infinite entanglement length. In 
fact, £e inherits the full critical behavior of £c, with 
the same critical exponent. All these results indicate 
again that multipartite entanglement plays a major role 
at quantum phase transitions, as argued in the previous 
paragraph. How important is the spontaneous symme- 
try breaking for this conclusion? In the next sections we 
develop this question for the 1-D XY spin chain. 



III. XY MODEL AND THE TWO-SPIN 
REDUCED DENSITY MATRIX 

The one-dimensional XY model in a transverse field is 
governed by the following Hamiltonian 

N N 
i=l i=l 

(3.1) 
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where erf , a = x,y, z, are the usual Pauli matrices. The 
model reduces to the quantum Ising model for 7 = 1 
and approaches the XX model [3{| as 7 — ► 0. The XX 
model belongs to a different universality class and we will 
therefore focus only on the parameter range < 7 < 1. 

The Hamiltonian l|3.ip is symmetric under a global ir 
rotation about the z axis (a x ^ — > —a x ^), which usu- 
ally implies a zero value for the magnetization in the x 
or y direction ((a x ^) = 0). However, as the magnetic 
field h is decreased (or J increased) this symmetry is 
spontaneously broken in the ground state (in the ther- 
modynamic limit) at A = J/h = Ai = 1, the first critical 
point (1CP). More specifically, the ground state is dou- 
bly degenerate with a finite magnetization ((<J X ) = ±M) 
in the x direction characterizing a ferromagnetic phase. 
It is also possible to define a symmetric ground state 
((a x ) = 0) using a superposition of the two degener- 
ate ones. Nonetheless, symmetric macroscopic states 
are just a theoretical construction with no physical exis- 
tence, since spontaneous symmetry-breaking mechanisms 
(superselection) rapidly destroy such coherent superpo- 
sitions (Schrodinger cats) in the thermodynamic limit 
[43 | . These unphysical states are called here symmet- 
ric states in contrast with the realistic broken-symmetry 
ones ((c x ) = ±M). Note that in the paramagnetic phase 
(A < 1) there exists no such distinction. 

By further decreasing the magnetic field a second phase 
transition occurs at A = A 2 (7) = — 7 2 , the second 

critical point (2CP). For magnetic fields smaller than this 
critical value, the correlation functions do not tend to 
their limiting value monotonically but in an oscillatory 
fashion [3^|. The Ising limit, 7=1, exhibits only the 
first critical point. 

As we will show in the following two Sections, for the 
calculation of the bipartite and multipartite (G(l) and 
G(2,n)) entanglement all we need is the reduced density 
matrix of two spins, a 4 x 4 matrix that can be expanded 
in tensor products of Pauli matrices and the identity er°: 



1 \ ^ a.B 



a,0 



where 



(3.2) 



(3.3) 



The reduced density matrix pij is obtained by tracing 
out all spins other than i and j. 

Remembering that pij is Hermitian with a unitary 
trace we are left with nine independent matrix elements 
for pi j , which are functions of the nine possible one and 



two-point correlation functions (pff = ',")■ This num 



0,a\ 



hJ 



ber can be further reduced by the symmetries of the prob- 
lem. In the XY model the global phase flip symmetry 
(global 7r rotation about the z axis) in the paramagnetic 
phase (A < 1) implies [ofcrj, pij] = 0, which imposes that 

(°i ) — ( a f a j) ~ ( a i a j) — 0) leaving only five inde- 
pendent correlation functions: (af},(a"a"),a = x,y,z, 



and (afOj). In the ferromagnetic phase (A > 1) this 
no longer holds since the Hamiltonian symmetry is not 
preserved by the ground state and we have to explic- 
itly evaluate the nine one and two-point correlation func- 
tions, (of), (o'i'o-'j), a = x,y, z and {o~^ v ^) were obtained 
in jH, Ha| . We are left then with three off-diagonal two- 
point correlation functions to calculate, (erf erf), (erf a?) 
and (of <rf). Finally, due to the translational symmetry 



of the model, pij depends only on the distance n = \i- 
between the spins, p"'- = p%>@ and p"'° = pf = p a is 
same for all spins. 
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We could be tempted to say that pij is real since the 



matrix elements of Hamiltonian l|3.ip are all real, and 
use this fact to eliminate p 3 ^ and pJJ* as both quantities 
appear in pij multiplied by the imaginary i. However, 
this argument can be misleading since this "symmetry" 
is not preserved in the ferromagnetic state in the ther- 
modynamic limit. As a counterexample consider for in- 
stance the Ising Hamiltonian with the nearest-neighbor 
coupling in the y direction and a magnetic field in the z 
direction (7 = —1). In the ferromagnetic phase we would 
have a finite value for (erf), which then results in complex 
elements in p^j. Fortunately, exact expressions for two 
of the three remaining off-diagonal correlation functions 
have been obtained by Johnson and McCoy [4fj. In fact, 
they have calculated the full time-dependent correlation 
functions (erf (0) erf (t)) and (erf (0)erf (t)). In particular, 
they have shown that at any time t, (erf (0)erf (t)) = 
for all values of 7 and h, which leads to p^ z = 0. For 
(erf (0)er?(t)), they have shown that the leading term for 
large n in both phases is linear in t, suggesting that p xy 
(t = 0) might be exactly zero. Numerical calculations of 
Pn V for small chains have confirmed that it does indeed 
vanish in both phases p9j . Collecting all the previous 
results the reduced two-spin density matrix pij can be 
written as 



/l + 2p z +plf p x +pf 



p- 
1 )J 



■Pij 
■Pi? 
-Pif 



1 



Pi : 
f p. 



Pi 



p* - p 



1.1 

xz 



p 

p-j 

1 



+Ptj 

: + p vv 

~ Pi* 



p- -p 



xz 

ij 



Pij 

p a 

p3 



P 



I J 

t xz 

ij 
xz 



1 - 2p z 



7 



+Ptj. 

(3.4) 

The last off-diagonal correlation function, p xz , was ob- 
tained in terms of cumbersome complex integrals in [13] 
rendering its explicit computation very tedious. How- 
ever, we were able to obtain bounds for it from the phys- 
ical restriction that all eigenvalues of pij must be posi- 
tive. Considering one of its eigenvalues as a function of 
results in a second-degree polynomial with negative 
second derivative (we have checked this for many values 
of A ranging from to 3 and for 7 ranging from 0.1 to 1). 
This allowed us to obtain tight lower and upper bounds 
for the value of p^ . This completes our construction of 
the reduced density matrix of two spins, which is all we 
need for the calculation of entanglement. In Figs. (H])- 
([5]) we plot the magnetization (erf) along a = x, z and 
the diagonal correlation functions, (0^0^), f3 = x,y,z, 
for nearest neighbors, j = i ±1, for later discussion. 
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Figure 1: (Color online) Magnetization along the a>axis for 
the XY model, with anisotropy. The first transition, 1CP, is 
apparent from the discontinuity at A = 1, for any anisotropy 
7- 




Figure 2: (Color online) Magnetization along the z-axis for 
the XY model. 



IV. BIPARTITE ENTANGLEMENT 

The bipartite entanglement between a pair of two- 
level systems (qubits or S = 1/2 spins) can be quan- 
tified using the concurrence C, since it is a monotonic 
function of the entanglement of formation 0| , a well es- 
tablished measure. The concurrence can be obtained 
from the density matrix of the two spins and is given 
by C = max{0, ei — e-2 — 63 — 64}, where ej, i = 1, . . .,4, 
are the square roots of the eigenvalues, in decreasing or- 
der, of the matrix R = pp. Here p = {a v ®a v )p*{a v ®a v ). 
Another pairwise measure of entanglement is the nega- 
tivity which is based on the Peres-Horodecki separability 




Figure 3: (Color online) Nearest neighbor (j = i±l) diagonal 
correlation function {afaf) for the XY model. 




Figure 4: (Color online) Nearest neighbor (j =i±l) diagonal 
correlation function (crfcrj) for the XY model. 

test [Hoi . IHlj . This test states that a separable state is 
always positive under partial transposition (PPT). This 
is also a sufficient condition for separability in the case 
of two-level systems. Thus, it is reasonable to quantify 
entanglement measuring "how much" the partially trans- 
posed density matrix is negative. A possible definition 
of negativity, which was proved to be an entanglement 
monotone [lfj], is given as N(n) = max{0, — 2min(Mfe)}, 
where Uk are the eigenvalues of the partial transpose of 
Pi,i+n and the label n denotes the distance between the 
qubits. The main advantage of the negativity over the 
concurrence is that the former is easier to compute than 
the latter. 

In the paramagnetic phase (A < 1) or for any value 
of A in the symmetric ground state (i. e., without sym- 
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Figure 5: (Color online) Nearest neighbor (j = i± 1) diagonal 
correlation function (ofcr?) for the XY model. 

metry breaking), the concurrence and the negativity are 
simple expressions in terms of the correlation functions. 
Mathematically, this is a consequence of the fact that 
the fourth-degree equations resulting from the diagonal- 
ization of R factorize into two equations of the second 
degree (H3| - The expression for the concurrence, valid 
for any system possessing the same symmetries as the 
symmetric ground state of the XY model is written as 

C{n) = max{0, C'(n), C"(n)}, (4.1) 

where 

C'(n) = 1(^-^1+^-1), (4.2) 

C"(n) = \ (\ P T+P y n V \ ~ V(l+P Z n Z ) 2 -HP Z ) 2 )^) 

The negativity expression derived exclusively for the XY 
model reads 

N(n) — max{0, -2 min[iii(n), u 3 (n)}}, (4.4) 

where 

ux{n) = -\ (l+pl* - vW+<) 2 +4(p z ) 2 )4.5) 

u 3 (n) = -l(l- p ™ +p yy-p* n *). (4.6) 

The expression for the concurrence could be written in 
this general form (independent of the particular values of 
the one and two-point correlation functions) because it 
is derived simply by imposing the positivity of pij and 
the fact that all eigenvalues of R are real numbers. We 
have also observed that for j 2 + 1/A 2 > 1, i.e. A < A 2 (7) 
it is always C'(n) and u 1 (n) that are relevant for the 
concurrence and the negativity through Eqs. (|4.ip and 



Figure 6: (Color online) Plot of the left-hand side of Eq. I|4.7p 
as a function of A: when the function is positive the con- 
currence for the XY model does not change upon symme- 
try breaking. The anisotropies are 7 = 1, 0.8, 0.6, 0.4, and 

0. 2 (from top to bottom). We have also plotted a small 
vertical bar to represent the position of the 2CP (A2 (7) = 

i/\/T^)- 

(|4.4p . respectively. On the other hand, for j 2 + 1/A 2 < 

1, i.e. A > A2 (7), it is C"{n) and u^{n) that appear 
in these two measures, respectively. The change from 
7 2 + 1/A 2 > 1 to 7 2 + 1/A 2 < 1 occurs at the 2CP, i.e. 

at A = A 2 (7) = Wl -7 2 - 

In the ferromagnetic phase (for the broken-symmetry 
state) the calculation of the concurrence is not so simple 
since we have a fourth-degree equation which factorizes 
into a first-degree one and a complicated third-degree 
equation. Although the latter can be solved exactly, the 
expressions for its roots are not very illuminating, render- 
ing a detailed general analysis unfeasible. Fortunately, it 
was demonstrated [H3| that for the Ising model the con- 
currence does not change upon spontaneous symmetry 
breaking. This opened the possibility for the use of the 
simple expression of the paramagnetic phase in the fer- 
romagnetic one. The analysis can be extended to the 
XY model since the reduced density matrices of the two 
models have a similar form. The condition for an iden- 
tical expression for the concurrence in the paramagnetic 
and ferromagnetic phases is 

V(l+P Z n Z ) 2 -±(P Z ) 2 + Pn ~ 2p V n V " 1 > 0. (4.7) 

In Fig. [6] we show the left-hand side of Eq. (|4.7j) as 
a function of A for the XY model. It can be seen that 
Eq. (|4.7j) always holds for the Ising model (7=1, first 
curve from top to bottom) but is violated after the 2CP 
(A > A 2 (7) = - 7 2 ) in the XY model (7 ^ 1). The 

critical value A2 (7) has been indicated by the vertical 
lines for each 7 in Fig. [6l 

Using the above expressions and the correlation func- 
tions depicted in Figs. <(T|) — CSj , including the bounds for 
p^ z obtained through the procedure explained in Sec- 
tion [nil we now analyze the bipartite entanglement be- 
tween any two spins of the XY chain. Notice that for the 
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Figure 7: (Color online) Lower bound of the concurrence for 
nearest neighbors obtained using the upper bound of pf z . In 
the limit of small magnetic field (large A) the entanglement 
decreases with increasing 7. Here 7=1, 0.8, 0.6, 0.4, 0.2, and 
0.1 from top to bottom in the A < 1 phase. 



Figure 8: (Color online) The lower and upper bounds for 
the nearest-neighbor concurrence (solid lines) and the con- 
currence in the symmetric state (dashed lines). Here 7 = 
1, 0.8, 0.6, 0.4, and 0.2 from top to bottom in the A < 1 phase. 
Note that for most of the anisotropics we can barely see the 
difference between the lower and upper bounds. 



following discussion we have calculated numerically the 
concurrence, not relying in the simplified formula l|4.ip . 
First, we have evaluated the concurrence of nearest neigh- 
bors (C(l)) for some values of the anisotropy 7. These re- 
sults are shown in Fig. where we plot the lower bounds 
for the concurrence. Moreover, in Fig. [8] we show both 
the lower and the upper bounds for the broken-symmetry 
state as well as the concurrence in the unphysical sym- 
metric ground state. As it is already known [6|,LZj], C(l) is 
not maximal at the 1CP. It is important to note that the 
bounds are very tight near the QPTs allowing us to cor- 
rectly characterize the behavior of the concurrence at the 
CPs. Indeed, Fig. [8] shows that the concurrence changes 
more abruptly (a diverging derivative) at the 2CP than 
at the 1CP where the spontaneous symmetry breaking 
occurs (A = Ai = 1)- We can also see that after the 1CP 
the concurrence starts to decrease, vanishing at the 2CP. 
This fact had already been observed (24| and it can be 
shown that at this point the ground state is completely 
separable [13] ■ 

Remarkably, the discrepancy between the symmetric 
and the broken-symmetry cases, in contrast to what one 
might expect, only occurs after the 2CP (A > A2), where 
the correlation functions tend to their limiting value in 
an oscillatory fashion [3S|. Thus, the spontaneous sym- 
metry breaking, which occurs already at the 1CP, has no 
influence on bipartite entanglement. Even after the 2CP 
the difference between the symmetric and the broken- 
symmetry states is small and becomes more pronounced 
only as 7 — -> 0, where the former has slightly less en- 
tanglement than the latter. Finally, the "origin" of the 
entanglement is different in the two states since in the 
symmetric case at the 2CP there is a change in the great- 
est eigenvalue of R, see Fig. [9] (this results in a change 
from C'(n) to C" (n) as the expression that contributes 
to the concurrence). This has been interpreted as a 
change in the kind of entanglement present in the ground 
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Figure 9: (Color online) Plot of the square root of the ei- 
genvalues of the matrix R (for 7 = 0.8), used to obtain the 
concurrence, for the symmetric case. We can observe a change 
in the greatest eigenvalue of R at the second critical point re- 
sulting in a change from C'(n) to C"(n) as the expression 
that contributes to the concurrence. 



state. However this is not true when one correctly em- 
ploys the broken-symmetry state since now it is the same 
eigenvalue of R that is maximal for all values of A, see 
Fig. [lOl We can also see that, as we approach the XX 
model (7 — ► 0), the concurrence decreases in the A < 1 
region, vanishing at the 1CP, and then increases in the 
ferromagnetic phase (A > 1). 

Similar conclusions can be derived for the negativity 
since the curves for the negativity as a function of A and 
7 are very close to the ones already shown for the con- 
currence. Here we only plot the negativity for nearest 
neighbors, N(l), in the symmetric state in Fig. [6l We 
should note that, as for the concurrence, the negativity 
of the symmetric ground state is close to the negativity of 
the broken-symmetry state in the ferromagnetic phase. 
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Figure 10: (Color online) Plot of the square root of the eigen- 
values of the matrix R (for 7 = 0.8), used to obtain the con- 
currence, for the broken-symmetry case. We can observe that 
for the broken-symmetry case there is no crossing of eigenval- 
ues of R at the second critical point. Actually, all eigenvalues 
vanish at the 2CP. This plot was obtained using the upper 
bound for p xz but there is no visible difference if we use the 
lower bound. 



N( 




Figure 11: (Color online) Nearest neighbor negativity in the 
symmetric ground state. We can see that at the critical point 
the entanglement decreases as we approach the XX model 
(7 - 0). 



We have also obtained the concurrence between next- 
nearest neighbors C(2) in the broken-symmetry state. In 
Fig. [32] we show its lower bound for some values of 7. We 
can see that in contrast to C(l), C(2) at first increases 
in the paramagnetic phase as we leave the Ising model 
in the direction of the XX model. As a function of A, 
C(2) reaches its maximal value just before A = 1 (the 
1CP), which increases as 7 — » 0. For small magnetic 
field (large A) we see the same behavior as for C(l): the 
entanglement increases as we approach the XX model 
(7 — > 0). The difference in the entanglement of the sym- 
metric and the broken-symmetry states, however, is more 
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Figure 12: (Color online) Lower bound for the concurrence 
for next-nearest neighbors using the upper bounds of the off- 
diagonal correlation function p% z ■ Here 7 = 1, 0.8, 0.6, and 0.4 
from bottom to top in the A < 1 phase. The broken-symmetry 
state was used in the ferromagnetic phase (A > 1). 
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Figure 13: (Color online) Concurrence for next-nearest neigh- 
bors in the broken-symmetry state at 7 = 0.2, obtained using 
the upper (triangles) and lower (squares) bounds of the off- 
diagonal correlation function p% z . For comparison, we also 
plot the concurrence in the symmetric state (crosses). 



pronounced now. In contrast to the broken-symmetry 
state, the entanglement in the symmetric state vanishes 
for A larger than a certain value. This is illustrated in 
Fig. [13] where we show the lower and upper bounds for 
the concurrence in the broken-symmetry state compared 
to the concurrence in the symmetric state for 7 = 0.2. 
We should note that for all values of 7 > 0.2, the bounds 
are tighter in comparison. 

The negativity for next-nearest neighbors has a similar 
behavior but with smaller values in the region of small 
fields (large A). This can be better viewed in Fig. fT4l 
where we plot N(2) in the symmetric state. Note that in 
the symmetric state the negativity (and also the concur- 
rence) vanishes for sufficiently small values of the field. 

We have also calculated the concurrence and the neg- 
ativity for spins three and four lattice sites apart, which 
show behaviors similar to C(2) and N(2). The only dif- 
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Figure 14: (Color online) Negativity for next-nearest neigh- 
bors in the symmetric ground state. We can see that at the 
first critical point (A = 1) the entanglement first increases as 
we leave the XX model (7 = 0) and then starts to decrease 
for sufficiently large 7. 



ference is that the entanglement is much smaller, decreas- 
ing as we increase the distance between the spins. It 
should be noted that the concurrence and negativity val- 
ues for next-nearest neighbors are significantly smaller 
than the ones for nearest neighbors, which means that 
bipartite entanglement is more concentrated on nearest- 
neighbor sites. 

One last fact about bipartite entanglement which we 
would like to mention is related to the origin of the non- 
analyticities of the concurrence and of the ground state 
energy per site £. It terms of the two spin reduced density 
matrix elements (pij) a f3 = (p)a/3 we have 



die 



and 



(7(1) = 2[(p) 4 i - (p) 22 ], 



where C(l) is the concurrence for the symmetric case and 
d\ stands for the second order derivative with respect to 
the tuning parameter A. As argued in [j| both d\£ and 
d\C(l) exhibit critical behavior through their dependece 
upon d\(p)22, since (p)u is well behaved. Looking at the 
correlation functions we observe that the divergence of 
d\{p)22 comes from the correlation in the z direction as 
(p)22 = 1 — p zz ■ However {p)a also has a dependence on 
p zz which cancels the one in (p)22, and we end up with 



d\e = --d xP z 



(4.10) 



Therefore, in terms of the correlation functions, the crit- 
ical behavior of df£ is originated on p z while the diver- 
gence of d\C(l) is given by p zz , since {p)a = p xx — p vv 
does not cancel the divergence of p zz . 



In sum, based on these measures of entanglement we 
can say that there is more bipartite entanglement around 
the 1CP when we approach the Ising model. When we 
approach the XX model the entanglement is more pro- 
nounced in the ferromagnetic phase. We expect this 
behavior to continue to hold for G(n) and N(n) when 
n > 4. It has already been pointed out that the range in 
which the concurrence has a finite value increases as 1 /j 
Q, being infinite for the XX model. Finally, we would 
like to stress that the concurrence does not change upon 
symmetry breaking in the Ising model. This fact shows 
that it suffers no influence during the symmetry-breaking 
process, although bipartite entanglement is able to mark 
the phase transition through a diverging derivative of the 
concurrence. 

In the XY model, on the other hand, the concurrence 
is different for symmetric and broken-symmetry states. 
However, this difference only appears at the second crit- 
ical point, not as a result of the symmetry breaking that 
occurs at the first critical point. Thus, we have shown 
that for both the Ising and the XY model, the symmetry- 
breaking phase transition does not affect the bipartite 
entanglement. In fact, as three of us have already ar- 
gued [42|, the fact that the concurrence does not depend 
on the magnetization in the x direction p x is a possible 
explanation for the fact that it is not maximal at the 
critical point. 



V. MULTIPARTITE ENTANGLEMENT 

After studying the entanglement between two spins in 
the chain we now focus our attention on the multipartite 
entanglement (ME) in the XY model as given by G(l) 
Q and G(2,n) |Il|, 52, H] . For a spin-1/2 chain these 
are given by 



(4, 



(4.9) and 



G(l) 



1 - 



1 N 



(5.1) 



G(2,n) 



d-1 



^ N-n 



+n> 



(5.2) 



where d is the dimension of the Hilbert space of the re- 
duced density matrices pj or Pj,j+ n , i.e., d = 2 for G(l) 
and d = 4 for G(2,n). For systems with translational 
invariance such as the XY model, these measures simply 
reduce to the the linear entropy of one spin and two spin 
n-sites apart, respectively. 

One advantage of the measures (|5.ip and (|5-2|) is that 
they are simple to evaluate since we just need the reduced 
one- and two-spin density matrices that were already ob- 
tained previously. This kind of measure has received the 
name of local entanglement measure / estimator since they 
depend on the reduced density operator of two spins, 
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which is a local quantity. Moreover, inspired by accumu- 
lated experience of many-body physics, we believe that a 
great deal can be learned from the knowledge of two-point 
correlation functions only. G(2, n) inherits the full non- 
analytical behavior of the elements of the reduced density 
matrix of two spins and it is possible to make a general 
link between divergences in the derivatives of the energy, 
which signal QPTs, and G(2,n) or its derivatives [44|. 
In other words, G(2,n) and/or its derivatives are able to 
signal QPTs as do bipartite entanglement measures, with 
the exceptions of the cases where the non-analyticities 
accidentally cancel out. In Ref. [Hj] it was shown that 
G(l) and G(2,?i) are maximal at the 1CP and zero at 
the 2CP, thus possessing the ability to map out the com- 
plete phase diagram of the XY model. Here we intend to 
analyze those results in vision of the differences between 
the symmetric and broken-symmetry cases/states. 

For the XY model Eqs. I|5.ip and (|5.2p above can be 
rewritten in terms of the one and two-point correlation 
functions, depicted in Figs. |1])-((I]), as 

G(l) = 1 - (a*) 2 - (a]) 2 (5.3) 

and 

G(2,n) = l-i(2(^) 2 + 2(a;) 2 + 2( ( 7^ + „) 2 + 

<^„> 2 + K^J 2 + K ( tj + „) 2 )(5.4) 

As we mentioned before, Eqs. l|5.3p and l|5.4p have already 
been shown to be maximal at the 1CP of the Ising [41| 
and XY [H[ models. We have argued that this behavior 
is due to the emergence of a finite value of the magnetiza- 
tion with the spontaneous symmetry breaking for A > 1 , 
since this is the correlation function in the expressions 
for G(l) and G(2,n) exhibiting the most abrupt change, 
despite being continuous at the 1CP. This will be made 
clearer as we investigate the behavior in the symmetric 
case. We also note that at the 1CP G(2,n) always in- 
creases as a function of n which is a strong indication 
of genuine ME. The von Neumann entropy (Entropy of 
Entanglement) of one spin has already been shown to be 
maximal at the 1CP @. 

We first show G(l) for the XY model as a function 
of A and 7 for the broken-symmetry case (see Fig. [T5)l . 
We can see that for all values of the anisotropy G(l) is 
maximal at the 1CP and decreases as one approaches the 
XX model. We have also checked that G(l) is zero at the 
2CP A = A 2 (7) = 1/yl — 7 2 as expected, since at this 
point the state is completely separable [Hj]. 

The upper bound for G(2, 1) for the broken-symmetry 
case is shown in Fig. [TU We note that it behaves sim- 
ilarly to G(l). To check the quality of the bounds we 
also plot the lower and upper bounds for G(2, 1) in the 
broken-symmetry case for three values of anisotropy in 
Fig. EH The same behavior was also found for G(2, 2) 
and G(2, 7) [II]. In fact, the value of G(2, n) for a fixed 
A always increases as a function of n, which is in con- 
trast to bipartite entanglement, and also in contrast to 




Figure 15: (Color online) G(l) for the XY model in the 
broken-symmetry case. 



G(2, 1 




Figure 16: (Color online) G(2,l) for the XY model in the 
broken-symmetry case. 



the correlation functions, which decrease as n increases. 
Therefore, G(2,n) must increase (see Eq. I5.4|) . either as 
a power law at the critical point or exponentially away 
from it. As mentioned before, this feature allows the def- 
inition of an entanglement length (see Ref. [H[ for more 
details) and indicates that at the CP the entanglement 
is more distributed in the chain than anywhere else. 

Since both G(l) and G(2,n) essentially show the same 
behavior, from Eq. I|5.3p we conclude that the magnetiza- 
tions (crj) and (erj) are the minimal quantities from which 
the multipartite entanglement over the chain can be in- 
ferred. Now we observe that near the XX model (7 — ► ) 
the ferromagnetic phase shows more entanglement than 
the CP, a feature observed for bipartite entanglement 
between two spins as well, as given by the concurrence 
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Figure 17: (Color online) Lower (solid lines) and upper 
(dashed lines) bounds for G(2,l) of the XY model in the 
broken-symmetry case for three values of the anisotropy: 
7 = 1, 0.6, 0.2. (from top to bottom in the A < 1 region). 



and negativity (see Figs. ([7]) and (JH|)). This behavior 
should be contrasted with Figs. (J])-©) f° r the magneti- 
zation and diagonal correlations. Notice that those one 
and two-point correlations (apart from (of £rj}) are in- 
variably smaller closer to the XX model than closer to 
the Ising model, in particular (crj) and showing the 
well know fact that, deep in the ferromagnetic phase the 
magnetization is destroyed by quantum fluctuations as 
ones goes from the Ising Hamiltonian to the XX limit. 
Thus it is reassuring to see that the proposed indica- 
tors of "quantum character", both bipartite and multi- 
partite entanglement measures, do indeed increase as we 
approach the XX model from the Ising one. This feature 
highlights the differences between classical and quantum 
correlations (entanglement) in the XY model. Thus, al- 
though the correlation functions involve both classical 
and quantum correlations, only a proper combination of 
them can reveal their entanglement content. 

We now compare the symmetric and broken-symmetry 
states. In Fig. [18] we have plotted G(l) and in Fig. [TjJ 
G(2, 1), both for two values of the anisotropy, 7 = 1 and 
7 = 0.4. It can be seen that in both cases the symmet- 
ric state (dashed line) does not show a maximum at the 
1CP. Instead, it is an increasing monotonic function of A. 
The crucial element here is that the magnetization in the 
x direction, (cr x ), and (erf a*) vanish in the region A > 1 
in the symmetric case, and the former is responsible for 
the difference in behavior between G(l) and G(2, 1). To 
check this, we have made (a x ) zero by hand in the expres- 
sion for G(2, 1) in the broken-symmetry case and verified 
that the result is very similar to the symmetric case. This 
indicates that the magnetization is the primary reason 
why G(l) and G(2, 1) are maximal at the 1CP and the 
results for the symmetric and broken-symmetry cases are 
different. The reader should remember that in the case 
of bipartite entanglement, given by the concurrence or 
the negativity, spontaneous symmetry breaking had no 
effect, only appearing to contribute in a different manner 
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Figure 18: (Color online) Comparison of G(l) in the symmet- 
ric (dashed lines) and broken-symmetry (solid lines) states 
for the XY model and for two values of anisotropy: 7=1 
(brown) and 0.4 (blue). 
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Figure 19: (Color online) Comparison of G(2,l) in the sym- 
metric (dashed lines) and broken-symmetry (solid lines) states 
for the XY model and for two values of anisotropy: 7=1 
(brown) and 0.4 (blue). 

after the 2CP. Furthermore, we observe that both G(l) 
and G(2,n) signal the two critical points in the broken- 
symmetry state, while the same is not true in the sym- 
metric state. The complete phase diagram could thus 
be drawn only by considering the non-analyticities of ei- 
ther of these two measures. Thus symmetry-breaking has 
more effect over multipartite entanglement. Again, the 
presence of the magnetization (ct x ) (which is highly sensi- 
tive to symmetry-breaking) in both G(l) and G(2, n) and 
its absence in the concurrence and negativity is respon- 
sible for the different behavior when one uses a broken 
symmetry ground state or not. 



VI. CONCLUSIONS 

We have extensively studied the entanglement prop- 
erties of the one dimensional XY model in a transverse 
magnetic field. We have in all cases assumed the chain 
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to be at T = and worked in the thermodynamic limit 
of an infinite chain. One of our goals was to charac- 
terize both the pairwise and multipartite entanglement 
of the ground state of the XY model. In order to do a 
complete analysis we were forced to consider two distinct 
ground states. The first one, which always preserves all 
the symmetries of the Hamiltonian, was called the sym- 
metric ground state. This state, however, is unphysical 
for A = JJh > 1 since in a realistic situation the global 
phase flip (global it rotation around the z axis) is always 
spontaneously broken. Therefore, we also considered a 
second state, namely the broken-symmetry ground state, 
where this symmetry no longer holds. 

For the broken-symmetry state we were able to show 
that, in contrast to pairwise entanglement, multipartite 
entanglement is maximal at the first critical point (where 
the XY model exhibits a diverging correlation length). 
This property is not observed in the symmetric state, 
in which case the multipartite entanglement increases 
monotonically as we decrease the external magnetic field. 
Furthermore, we have also shown that the concurrence 
does not change in the vicinity of the symmetry-breaking 
critical point whether we employ the symmetric or the 
broken-symmetry state. On the other hand, we have 
explicitly shown that the multipartite entanglement de- 
pends strongly on the symmetry of the ground state. 
This result suggests that, as is the case for the XY model, 
the behavior of the multipartite entanglement may be 
intimately connected with the spontaneous symmetry 
breaking mechanism. We should also remark that only 
after the second critical point is the concurrence depen- 
dent on which state we use, and this is probably because 
the symmetry-breaking has a more pronounced effect in 
the oscillatory behavior the correlations show after this 
point. 

We have arrived at another interesting result by notic- 
ing two important facts. First, for spins three or more 
sites apart there exists no pairwise entanglement what- 
soever @, 0|. Second, we have shown that multipartite 
entanglement is never zero at the first critical point (be- 
ing maximal for the broken-symmetry state) . Combining 
these two facts we are led to conclude that the long range 



correlations at this critical point are a consequence of the 
existence of multipartite entanglement and not pairwise 
entanglement. Moreover both multipartite and bipartite 
entanglement tend to be larger than at the 1CP in the 
ferromagnetic phase as one approaches the XX model 
(7 — ► 0), increasing monotonically after the 2CP. This 
feature contrasts with the one and two-point correlations 
that tend to decrease as 7 — > 0, showing that the en- 
hanced quantum correlations (entanglement) in this re- 
gion can only be appreciated through a proper combi- 
nation of the correlation functions. Since the XX model 
belongs to a different universality class, nothing can be 
said about entanglement at 7 = from our study. It 
would be certainly interesting to investigate how entan- 
glement develops in the distinct phases of the XX model. 

Finally, our results have shown that the entanglement 
contents of the two possible ground states, i.e. the 
symmetric and the broken-symmetry state, are different. 
Therefore, it is of utmost importance to clearly and ex- 
plicitly verify which state one is using in any entangle- 
ment analysis made for the XY and related models in 
order to avoid any possible confusion. 

Note: Upon finishing this manuscript we became aware 
of a work where the effects of s ymm etry breaking on the 
concurrence are also addressed [53| , and a detailed study 
of the concurrence and the entanglement between one site 
and the rest of the chain (one-tangle) • 
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